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Abstract 

The existence probability Ep and the percolation probability P of the bond percolation on 
rectangular domains with different aspect ratios R are studied via the mapping functions between 
systems with different aspect ratios. The superscaling behavior of Ep and P for such systems with 
exponents a and 6, respectively, found by Watanabe, Yukawa, Ito, and Hu in [Phys. Rev. Lett. 93, 
190601 (2004)] can be understood from the lower order approximation of the mapping functions 
fu and QR for Ep and P, respectively; the exponents a and b can be obtained from numerically 
determined mapping functions fa and qr, respectively. 
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I. INTRODUCTION 



Universality and scaling in critical systems have attracted much attention in recent 
decades It has been found that critical systems can be classified into differ- 

ent universality classes so that the systems in the same class have the same set of critical 
exponents According to the theory of finite-size scaling 0, 3, 4|, if the dependence of a 
physical quantity M of a thermodynamic system on a parameter t, which vanishes at the 
critical point t = 0, is of the form M{t) ~ near the critical point, then for a finite system 
of linear dimension L, the corresponding quantity M{L, t) is of the form 



M{L,t) 



-ayt 



F{tL 



(1) 



where yt (= z/^^) is the thermal scaling power, F{x) is the scaling function, and x = tV^^ 
is the scaling variable. It follows from ([T]) that the scaled data M{L,t)L°'y^ for different 
values of L and t are described by a single finite-size scaling function (FSSF) F{x). Thus 
it is important to know general features of the scaling function under various conditions. 
On the basis of the renormalization group arguments, Privman and Fisher j^, 4] proposed 
that systems in the same universality classes can have universal finite-size scaling functions 
(UFSSFs). 



Percolation models 



Id, 



are ideal systems for studying universal 



and scaling behaviors near the critical point. The system is defined to be percolated when 
there is at least one cluster extending from the top to the bottom of the system. The average 
fraction of lattice sites in the percolating clusters is called the percolation probability [s], [g] 
(also known as percolation strength) and will be denoted by P in the present paper. The 
probability that at least one percolating cluster exists in the system is called the existence 
probability jol or crossing probability [isl, and will be denoted by Ep in the present paper 



Using a histogram Monte Carlo method and other Monte Carlo methods, Hu and 
collaborators found that the FSSFs depend sensitively on the boundary conditions [l^ 
and shapes 1^ of the lattices and many two-dimensional percolation models (including 
bond and site percolation models on square, honeycomb, and triangular lattices, continuum 
percolation of soft disk and hard disks, bond percolation on random lattices, etc) [ij] can 
have universal FSSFs (UFSSFs) for their Ep and P, and many three-dimensional percolation 
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models 
lattice) 



including site percolation on sc, bcc, and fee lattices and bond percolation on sc 



18l | have universal critical exponents and UFSSFs for their Ep and P. 



In the present paper, we will study the bond percolation model on Li x L2 rectangular 



lattices (domains) 19|, |20|, especially the dependence of Ep and P on the aspect ratio R, 
which is the ratio of the horizontal length Li to the vertical length L2 = L of the system, i.e., 
R = L1/L2. For this system, it is well known that the critical bond density pc (also called 
bond occupation probability) is 1/2, the critical exponent (3 for the percolation probability 
P is 5/36, ?/t=l/z/ = 3/4, and the critical exponent for the existence probability Ep is 
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2l| , which and Eq. ([T]) imply that the dependence of Ep on the the linear dimension L, 



the bond density p, and the aspect ratio R can be written as 

Ep{L,p,R) = FR{eLy^), [e ^ \p - p,\l p,) (2) 

near the critical point Or jsl, [ll| . The function Fr is universal for various percolation models 
considered in Ref. [l7|. The aspect ratio dependence of the function Fr is not clear; e.g., 
the systems with close values of R should have similar forms of Fr, however, Eq. ([2]) do 
not tell us how similar they are. While the value of the existence probability at the critical 
point Ep{pc) was obtained by Cardy [20], such critical Ep approaches quickly to 1 as the 
aspect ratio is increased, e.g., it is difficult to distinguish the value for R = 16 from that for 
i? = 32. Additionally, the region for the rapid increase of Ep shifts to the smaller value of 
p with larger R and fixed L, which cannot be described quantitatively. Therefore, we need 
other kinds of approach to study the aspect ratio dependence of the UFSSFs. 

Recently, we found "superscaling" of the existence probability Ep and the percolation 



probability P for bond percolation on rectangular domains with different aspect ratios 22] : 
Ep{L, p, R) ~ F{e'Ly'R^) and P(L, p, R) ~ {Ly'R'')-^F{e'Ly'R^), with new exponents a and 
b, where e' = {p — p^) with p'^ being the effective critical point; a and b were determined by 
a fitting procedure to be a = 0.14(1) and b = 0.05(1), respectively. In the present paper, 
we define mapping functions to connect UFSSFs for percolation on domains with different 
aspect ratios and find that a and b can be obtained from numerically determined mapping 



functions fn and gR for Ep and P, respectively [23 1. 
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II. MAPPING FUNCTION AND FINITE-SIZE SCALING 



In order to make the meaning of the mapping function becomes clear, we first introduce 
the mapping function for the percolation model on the square lattices with the same aspect 
ratio. Consider two different square systems A and B with the linear sizes and Lb{> La), 
respectively. The existence probability of the system B at density (also called occupation 
probability) p is denoted by Ep{LB,p)- The existence probability is a monotonic increasing 
function whose value changes from to 1 as the density increases from to 1. Therefore, 
for every Ep{LB,p) of system B, there is a corresponding value p' of system A so that the 
existence probability of system A at p', denoted as Ep{LA,p'), is equal to Ep{LB,p) and we 
have 

Ep{LB,p) = Ep{LA,p')=EpiLAjQ{p)), (3) 
which defines the renormalization group transformation from p to p' in Ref . [9|] . The mapping 
from p to p' is represented as the mapping function 

P = fA^Bip) ^ fgip). (4) 

The mapping function depends only on the ratio of the system size Q = Lb /La, and 
therefore, we denote the mapping function by fg. The critical point of the system, p^, can 
be determined from the fixed point of Eq. ([3]) : 

EpiLB,pc) = Ep{,LA,Pc), 

which is equivalent to 

/q(Pc) = pc- (5) 
For fixed Q, Pc of Eq. ([5]) approaches the exact critical point 1/2 for the bond percolation 
on the square lattice as La — > oo. Table 1 of the first reference in shows that pc of Eq. 
([5]) is already very accurate for = 16 — 20 and La = 8 — 16. 

The expansion of the mapping function at the critical point can be written as 

P = /q(p) (6) 
dfQ 



/q(Pc) + 



dp 



(P-Pc) + ---. (7) 



Pc 

It is well known that the thermal scaling power yt can be obtained from the equation 



which imphes that the differential coefficient of the first order in Eq. ([7]) is given by 



Q^*. (9) 



dp 

We can also have the above relation from the viewpoint of the mapping |24l ]. 
Using Eqs. ([7]), ([5]) and (111), we can rewrite Eq. ([3]) as 

Ep{LB,p) = Ep{LA,Pc + Q'''{p- Pc)) 

= Ep{LA,Pc + ^L~y') (10) 

with the scaling variable e = (p — Pc)L^B- Since the system sizes La and Lb are arbitrary, 
we can choose La as unity and rewrite Lb a.s L. Then we have 

E,mp) = E.p{l,p, + e) = F{e), (11) 

with the scaling function F. Equation (ITTl) implies that the finite-size scaling theory corre- 
sponds to the first order approximation of the expansion in Eq. ([7]). 

To see the behavior of the mapping function fq, we determine /q from numerically 
obtained existence probabilities. The procedures to determine fq is as follows, (i) Prepare 
graphs of the existence probabilities vs. density, (ii) Draw a line parallel to the horizontal 
axis, (iii) Determine the two intersection point p and p' (see the inset of Fig. [T](a)). (iv) Plot 
all the pairs (p, p') by sweeping from Ep = to 1, then the function p' = /q(p) is obtained. 
The calculated mapping functions of the bond percolation on square lattices are shown in 
Fig. [H^a) with aspect ratio R = 1 and system sizes L = 64, 128 and 256. Free boundary 
conditions are taken for this and other systems in the present paper. All of these functions 
have linear forms and have the single intersection point, namely the critical point. The ratio 
of the slope of these lines gives the value of the critical exponent u = 1/yt. 

Similar arguments can be applied for the percolation probability P. A mapping function 
gg is defined by, 

L^^y^P{LB,p,R) = L^/^P{LA,gQ{p),R) (12) 

with a system size ratio Q = Lb /La- Here, we defined the mapping function with the 
additional factors L^^ and L^^ on the left-hand and right-hand sides of Eq. (fT2l) . With 
such factors, both sides of Eq. f|T2l) are proportional to the FSSF for P as can be seen from 
Eq. dl]). For Ep, the critical exponent is |21l] and we need not add such factors in Eq. ([3]) 
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to define the mapping function. Tlie functions qq are shown in Fig. [T](b). Figure 1 shows 
that two kinds of functions /g and qq are identical near the critical point. 

The UFSSFs corresponding to mapping functions in Fig. [T] are shown in Fig. [2J The 
goodness of the scahng in these figures come from the fact that the mapping functions are 
well approximated by the linear function when the systems have the identical aspect ratio. 
For the systems with different aspect ratios, the mapping functions for Ep becomes highly 
non- linear function as shown in Fig. 3(b) below. Therefore, the scaling with low-order 
approximation does not show good data collapse for large values of the scahng variables (see 
Fig. 4 in Ref. [2^). 



III. MAPPING FUNCTION AND SUPERSCALING 



As shown in Fig. 3(a), we define the mapping function //j(p) to map the existence 
probability of the bond percolation model on RL x L lattice, Ep{L, p, R), into the existence 
probability of the bond percolation model on L x L lattice, Ep{L, fni^p), 1), with following 
equation 

Ep{L,p,R)=Ep{L,fR{p),l). (13) 

The curves of the existence probabilities for the same L and different R do not intersect, 
since their values at the point pc of the thermodynamic system are different j2o| . Therefore, 
we introduce the effective critical point p[, which satisfies the equation 



Pc 



and consider the expansion of the mapping function at p[ as 



fRip) = mp'c) + E 



R 



^ dp" 

71=1 



{p-p'cT 



Pc + 



dp 



(P - p'c) + 



dp2 



(P - p'c) 



+ 



(14) 
(15) 



The expansion coefficients depend only on R. With similar arguments in Ref. ^] , the aspect 
ratio dependence of the differential coefficient of the first order can be assumed to be 

d/ij 



dp 



R\ 



(16) 
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and taking into account up to the first order derivative of tlie expansion, we obtain tlie 
approximated mapping function as, /i?(p) — pc oc W'i^p — p^). From tfie finite-size scaling, 

E,m Up), 1) ~ f iifnip) - Pc) Ly^) , (17) 

and Eq. (fT3|) . we can derive tlie superscaling form to be, 

E,mp,R) = E,{LJM,1) (18) 

~ Fiifnip)-Pc)Ly^) (19) 
~ Fiip-p'jLy^R'^). (20) 

The mapping functions fn of existence probabilities of the bond percolation model on 
the square lattices are shown in Fig. [3](b) for the system size L = 256 and aspect ratios 
R = 1,2,4,8 and 16. It shows that these functions do not have any intersection points, 
since Ep{L,p,R) > Ep{L,p,l). Additionally, while the mapping functions of fg in Fig. [1] 
have linear forms, the curves of in Fig. 3(b) for R> 1 are not linear. 

Consider the mapping of the percolation probabilities with different aspect ratio as, 

P{L,p,R) = P{L,gn{p),l), (21) 

with a mapping function g^i. Similar to the case of the existence probability, the behavior 
of the differential coefficient of qr can be assumed as 



dp 



R\ (22) 

P'c 



with the effective critical point p'^ defined by Qnip'^) = Pc as shown in Fig. 4(a). Note that, 
the value of effective critical point is different from the one defined by /i?(Pc) = Pc- We can 
obtain the superscaling formula as, 

P(L, p, R) ~ L'^y^F ((p - p',)Ly^R'') . (23) 

This scaling form is different from the form 

P(L, p, R) ~ {Ly^R'r^F ((p - p;)L^'i?^) (24) 



obtained by a heuristic argument in 22| . The mapping functions of percolation probabilities 



are determined and shown in Fig. HJ^b). While they also do not have any intersection points, 



they have almost hnear forms. The scahng plot based on Eq. (123|) is shown in Fig. [51 It 
shows good scaling behavior especially around at the effective critical point. Please note 
that the right-hand sides of Eqs. fl23p and fl24p differ only by the factor R~^^ which is very 
close to 1. 

Since the mapping functions //j for Ep have finite curvature, the scaling formula 
Ep{L, p, R) ~ F{e'Ly^R°') will be improved if we take into account higher order. This 



nonlinear scaling form was shown in [2^, and it corresponds to approximate mapping func- 
tion with a quadratic form as /h(p) = C2p^ + Cip + Cq, which is the expression obtained by 
taking into account up to the second order derivatives in Eq. ( |T5l) . With the coefficients 



C2,Ci and Cq listed in [22|], we calculate the differential coefficients and show the results in 
Fig. ini^a). The mapping functions of the percolation probabilities have almost linear form, 
we just fit them with linear function to obtain the differentiations. The obtained differential 
coefficients of the percolation probabilities are shown in Fig.[6](b). Figures [6](a) and (b) show 
power law behavior i?*^ and R!^ with values of exponents a = 0.14 and h = 0.05 which are 



consistent with the results in 22]. Thus a and b can be obtained from mapping functions 



fR and gji for Ep and P, respectively. 

IV. SUMMARY AND DISCUSSION 

In summary, we have studied the aspect ratio dependence of FSSFs of existence and 
percolation probabilities based on the idea of the mapping functions. From the analysis 
on the numerically obtained mapping functions, the superscaling behaviors are found to be 
the first order approximation of the mapping functions at the effective critical point. Since 
the mapping regime used here is pretty general, we can apply it for any other systems and 
geometric conditions. 

Many lattice phase transition models, such as the Ising model, the Potts model, the dilute 
Potts model, the hydrogen bonding model for water, lattice hard-core particle models, etc. 



have been shown to be corresponding to some correlated percolation models 25|]. Some of 



such correlated percolation models have been found to have good finite-size scaling behaviors 



261]. It is of interest to use mapping function^ to study super-scaling in such correlated 
percolation models and other critical systems j27 ]. 

Another interesting development for critical systems is finite-size corrections for lattice 
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Ising model on finite lattices 



phase transition models [28|, l29l. l30l . l3ll . |32| . |33| . Using exact partition functions of the 



3J] and exact finite-size corrections for the free energy, the 



internal energy, and the specific heat of the Ising model 32], Wu, Hu and Izmailian 35] 
obtained UFSSFs for the free energy, the internal energy, and the specific heat with analytic 
equations, which are free of simulation errors. Based on finite-size corrections for the dimer 
model on finite lattices [36|, Izmailian, Oganesyan, Wu, and Hu (37| obtained UFSSFs for 



the dimer model with analytic equations, which are also free of simulation errors. It is of 
interest to extend above results for the Ising and dimer models to rectangular domains with 
various aspect ratios. We can then study superscaling behavior of physical quantities of the 
Ising and dimer models with the help of mapping functions. 
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FIG. 1: Mapping functions of the bond percolation model on square lattices with free boundary 
conditions (a) fg for the existence probabilities Ep as & function of p, and (b) gg for the percolation 
probabilities P as a function of p. The insets in (a) and (b) show how to construct the mapping 
functions for Ep and P, respectively. 



12 




1.2 
1 

0.8 
0.6 
0.4 
0.2 




1 1 

(b) 


1 1 


1 1 




/ Q= 


=1 — 




/ Q= 


=2 


— 1 1 


1 1 


=4 

1 1 



-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



FIG. 2: UFSSFs of (a) existence probabilities and (b) percolation probabilities corresponding to 
Fig. [TJ The aspect ratio is -R = 1, and the system sizes are L = 64((5 = 1), 128(Q = 2), and 
256((5 = 4), respectively. The scaling variable is denoted by e = (p — pc)LyK 
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FIG. 3: (a) Procedures to obtain the mapping functions of existence probabilities with different 
aspect ratios, (b) Obtained mapping functions of systems with L = 256 and = 1,2,4, 8 and 16 
(from bottom to top) . The sohd circles are the effective critical points p'^,. 
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FIG. 4: (a) Procedures to obtain the mapping functions of percolation probabilities with different 
aspect ratios, (b) Obtained mapping functions of systems with L = 128 and i? = 1,2,4, 8 and 16 
(from bottom to top). 
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FIG. 6: Plots of (a) u = d//dp|p/_ and (b) v = dg/dp\pi^ as a function R. The solid lines are i?" 
and Ef' with exponents a = 0.14 and b = 0.05, respectively. Decimal logarithms are taken for both 
axes. 
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